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ABSTRACT 

The ALHAMBRA survey aims to cover 4 square degrees using a system of 
20 contiguous, equal width, medium-band filters spanning the range 3500 A to 
9700 A plus the standard JHKs filters. Here we analyze deep near-IR number 
counts of one of our fields (ALH08) for which we have a relatively large area 
(0.5 square degrees) and faint photometry (J=22.4, H=21.3 and K=20.0 at the 
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50% of recovery efficiency for point-like sources). We find that the logarithmic 
gradient of the galaxy counts undergoes a distinct change to a flatter slope in 
each band: from 0.44 at [17.0, 18.5] to 0.34 at [19.5, 22.0] for the J band; for 
the H band 0.46 at [15.5,18.0] to 0.36 at [19.0,21.0], and in Ks the change is 
from 0.53 in the range [15.0,17.0] to 0.33 in the interval [18.0,20.0]. These 
observations together with faint optical counts are used to constrain models that 
include density and luminosity evolution of the local type-dependent luminosity 
functions. Our models imply a decline in the space density of evolved early-type 
galaxies with increasing redshift, such that only 30% - 50% of the bulk of the 
present day red-ellipticals was already in place at z ~ 1. 

Subject headings: cosmology: observations — galaxies: photometry — galaxies: 
evolution — surveys — galaxies: high-redshift — infrared: galaxies 



Introduction 



It is well understood that the stellar masses of galaxies are better examined with near- 
IR (NIR) observations compared to shorter wavelengths mainly because the near-IR light 
is relatively less affected by recent episodes of star formation and by internal dust extinc- 
tion. Moreover, the K-corrections are also smaller and better constrained in the NIR and 
hence massive high redshift objects are relatively prominent in the NIR. Despite this rel- 
ative insensitivity to luminosity evolution and the e ffects of dust, the hope of using the 



NIR counts to constrain the cosmological para meters ([Gardner et al.lll993l ; iMoustakas et al. 



19971 : iHuang et al.l l200ll ; iMetcalfe et al.l 120011 ) has not proved feasible because evolution in 
the space density of galaxies was soon un derstood to be of compa rable significance for the 
NIR counts as the cosmological curvature (IBroadhurst et al.lll992h . 



Disentangling the effects of cosmology from evolution is not straightforward even in the 
NIR, and now it has become more appropriate to turn the question aro und and make use o f 
the impressive constraint s on the cosmological parameters from W MAP (jSpergel et al.ll2003l ). 



and type la supernovae (IRiess et al 



carefully the rate of evolution (jMartini 



1998 



Per 



2001 



mutter et al.lll999l). in order to measure more 



Saracco et al. 12001: Cristobal-Hornillos et al. 



20031 : lEliche- Moral et al.ll2006l ). In addition, imaging in the NIR has progressed well with 



1 Based on observations collected at the German-Spanish Astronomical Center, Calar Alto, jointly oper- 
ated by the Max-Planck-Institut fur Astronomie Heidelberg and the Institute de Astroffsica de Andalucfa 
(CSIC). 
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fully cryogenic wide-field imagers now available on several 4m class telescopes. We also have 
at our disposal now much better estimates of the luminosity functions of different classes of 



galaxies from the larg e local redshi ft surveys in par ticular the SDSS (IBlanton et al 
Nakamura etUboO j ). or 2dFGRS JCole et al.ll200lh . 



2003 



The evolution of the luminosity functions has been addressed making use of spectro- 
scopic redshift surveys. However, the results of those studies differ due to the limitations 
in terms of low st atistic, or small fields probed which lead to uncertainties due to the large 
scale structure. In lFontana et al.l (12004 ) a mild (20-30%) evolution in the number density of 
massive objects since z ~ 1 is found. A lso a roughly constant number density for red galaxies 
to z = 0.8 is foun d in llm et al.l (120021 ). Using COMBO-17 photometric re dshift information 
(I Wolf et al.ll2003l ). and the rest-frame color bimodality at each redshift, iBell et al.l (120041 ) 
with a sample of ~25000 galaxies, concluded that the stellar mass in red-galaxies has in- 
crease in a fa ctor 2-3 from z ~ 1 to the present. Combining spectroscopic with photometric 
redshift data lllbert et al.l (120061 ) point to an increase in a fac tor of ~ 2.7 for th e density 
of red bulge-dominated galaxies between z = 1 and z = 0.6. iFaber et al.l (120071 ) found a 



different evolution since z ~ 1 in number densities in the red and blue galaxy population, 
being constant for the blue galaxies, while the number density of the red galaxies increases 
in a factor 3. Wide field imaging with larger covered area, and greater numbers selected to 
uniform faint limits is complementary to the redshift surveys in examining statistical models 
proposed for evolution. 

In practical terms it is most useful to combine faint NIR counts with deep blue counts 
when examining models of evolution to contrast the effects of lu minosity and dens i ty evo - 
lution which af f ect th ese two spectral ranges in different ways. In iMcCracken et al.l (120001 ); 
Metcalfe et al.l (120011 ) the authors use non-evolving models with a higher 0* normalization 
in the B-band even if this leads to an over-prediction of bright galaxies than is observed. 
Huang et al.l (120011 ) pointed out that both the optical and NIR counts present an excess 
over the no-evolution models, finding passive evolution models more suitable to match the 
distributions. The authors emphasize nevertheless their disappointment with the fact that in 
the passive evolution models the faint number counts are dominated by early-type galaxies, 
whereas the real data show that spiral and Sd/Irr galaxies are the main contributors to the 
faint counts even in the K-band. 



A characteristic feature of the NIR galaxy counts reported in several wor ks (I Gardner et al 



19931 ; iHuang et al.l 1200 ll ; ICristobal-Hornillos et al.ll2003l ; llovino et al.l 120051 ) is the change of 
slope at 17 < Ks < 18. This distinctive flattening is not observed in the B-band counts. 
This effect has been interpreted in terms of a change in the dominan t galaxy population , 



becoming increasingly dominated by an intrinsically bluer population ([Gardner et al.l 11993 
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Eliche-Moral et al.ll2006l ). In the model described in lCristobal-Hornillos et al.l (120031 ) a delay 
in the formation of the bulk of the early-type galaxies to Zf orm < 2 and the presence of a 
dwarf star-forming population are invoked to match the Ks-band counts. A similar dwarf 
st ar-forming population at z > 1, that is not present at lower redshift was found compatible 



in 



Metcalfe et al.l (119961 . 1200 ll ) but that work uses a go = 0.5, A = 0.0 cosmology, requiring 



some revision. 



Alternatively, an increase of 0* for late-type galaxies, driven via mergers, can produce 
similar results without in troducing an ad hoc population that is unseen in the local LFs 
(jEliche-Moral et al.ll2006l ). In any case, a low Zf orm ~ 1.5 for the ellipticals remains necessary 
to generate a significant decrease in the number of red galaxies and to account for the 
distinctive break in the NIR count slope at Ks~17.5. 

Here we use the NIR data from t he first comp l eted A LHAMBRA field, hereafter ALH08 
(details of the project can be found in lMoles et al.l (120081 ) and http://www.iaa.es/alhambra]). 
The limiting magnitudes (at S/N=5 in an aperture diameter of 2xFWHM) reached in the 
3 NIR bands are in mean for the eight frames in ALH8: J=22.6, H=21.5 and Ks=20.1 with 
a 0.3 rms (Vega system), and the total area covered amounts to ~ 0.5 square degrees. The 
completed survey will cover 8 independent fields with a total area of 4 square degrees. The 
ALHAMBRA survey occupies a middle ground in terms of the product of depth and area 
in all three standard NIR filters. The bright end of the counts is well constrained by our 
relatively large area allowing a careful examination of the location and size of the break in 
the count-slopes in J,H and Ks at fainter magnitudes. We have paid special attention to 
S/G separation, which at the intermediate magnitude range, is effectively achieved using 
optical-NIR color indices by combining our ALH08 NIR data with the corresponding Sloan 
DR5 data. 

Unless specified otherwise, all magnitudes here are presented in the Vega system, and 
the favored cosmological model, with H = 70 kms _1 Mpc _1 , Qm = 0.3, Qa = 0.7 was 
adopted through this paper. 



2. Observing Strategy and Data Reduction 



The ALHAMBRA survey is collecting data in 23 optical-NIR filters using the Calar 
Alto 3.5m telescope with the cameras LAICA for the 20 optica l filters, and OMEGA2000 
for the NIR, JHKs filters (IMoles et al.l 120081 : iBenitez et al.l 120081 ). In this paper we discuss 
the galaxy number counts in the J, H and Ks bands computed in the completed ALH08 
pointing. Due to the OMEGA2000 and LAICA geometries, two parallel strips of ~ 1 degree 
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x0.25 degrees are acquired in each of the eight ALHAMBRA fields. Each of these strips 
is covered by four OMEGA2000 pointings. Fig. [1] shows the central part of the processed 
image for one such pointing in ALH08. 

The OMEGA2000 camera has a focal plane array of type HAWAII-2 by Rockwell with 
2048 x 2048 pixels. The plate scale is 0.45 arcsec/pixel, giving a full field of view of ~236 
arcmin 2 . The JHKs images were taken using a dither pattern of 20 positions, with single 
images time exposures of 80 s in J, 60 s in H and 46 s in Ks band (obtained using respectively 
16, 20 and 23 software co-adds). The total exposure time was 5 ks in each of the three filters. 

Due to the high quantity of raw data images we have implemented a dedicated reduction 
pipeline to process the images. This code will be presented and discussed in a forthcoming 
technical paper. The use of the reduction pipeline guaranties the homogeneity of the process, 
and allows us to perform a first automatic analysis of the resulting images and to verify 
the quality of the products in real time. The magnitude at the 50% of recovery efficiency 
for point-like sources in the eight final images corresponds in average to J= 22.4 ± 0.24, 
H= 21.3 ± 0.14 and Ks= 20.0 ± 0.13 in the Vega system. The galaxy number counts have 
been computed in the high signal to noise region of the final images, covering a total of 
~1600 arcmin 2 , or 0.45 sq deg. 



2.1. Flat-fielding and sky subtraction 



Firstly, the individual images of each observing run were dark-subtracted, and flat- 
fielded by super-flat images constructed with the science images in each filter. In the case of 
NIR imaging it is specially important to remove properly the high sky level that is changing in 
short timescales. The sky structure of each image was removed with the XDIMSUM package 
( [Stanford et al.lll995l ) with a sky image constructed with the median of the 6 images closest 
in acquisition time, which in the case of the J, H, Ks filters correspond to timescales of 480, 
360, 276 s respectively. During this process cosmic ray masks for each individual image were 
also created. 



SExtractor (IBertin fc Arnoutslll996l ) was used with the preliminary sky subtracted im- 
ages to compute the number of detected objects (above a given S/N and with ellipticity 
lower than a certain limit), and to make a robust estimate of the FWHM of each image. 
This information together with the sky level variation was used to automatically remove low 
S/N images, lying outside the survey requirements, and also to identify exposures presenting 
problems such as telescope trailing. 
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Fig. 1. — Central part (~ 9.2' x 6.4') of one of the pointings in the ALH08 field in the Ks 
band. 
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2.2. Ghost and linear pattern masks 

The readout of the detector array produces ghost images coming from bright stars. As 
those spurious effects are replicated in all the readout channels (separated by 128 pixels) of 
the detector's quadrant where the real bright source is located, it was possible to mask them 
before image combination. 

Linear patterns produced by moving objects were located in the images using two dif- 
ferent approaches: i) Objects with high ellipticity and pixel area were identified from the 
individual image SExtractor catalogs, ii) Linear patterns split in multiple spots and struc- 
tures were located using the Hough transform. Linear patterns detected by any of the two 
methods were masked. 



2.3. Image combination 

Using XDIMSUM, the images were combined masking the cosmic rays, bad pixels and 
linear patterns. Those preliminary co-added images were used to create object and bright 
cores masks. Object masks were used to cover the sources in the second iteration of the 
XDIMSUM sky subtraction procedure. Bright cores masks operate when constructing an 
improved version of the cosmic ray masks. 

At this point, the individual images have been dark current corrected, flat fielded and 
sky subtracted. Each individual image has an associated mask containing the bad pixels, 
cosmic rays, linear patterns and ghost images. 



To perform the final combination of the processed images we used SWARP (IBertin et al. 



20021 ) . a code that combines the images, correcting at the same time the geometrical distor- 
tions in the individual images using the information stored in WCS headers. The astrometric 
calibration of each individual image and the update of its WCS headers was done by the 
pipeline using an automatic module. During the SWARP combination the extinction varia- 
tions among the images were also corrected. More technical details on the image combination 
are given in appendix [A] 



2.4. Photometric calibration 



For the photometric calibration we used the 2MASS catalogs (ICutri et al.ll2003l ). After 
having combined all the images of each pointing, the objects in common with the 2MASS 
catalogs were found and those with higher signal to noise ratio selected to compute a zero 
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point offset. 

In Fig. [2] an example of the photometric calibration for one of the pointings in ALH08 is 
shown. The histogram with the rms in the photometric calibration for the ALH08 pointings 
in the three NIR bands is shown in Fig. [31 The mean value for the rms is 0.028 ±0.006 mags, 
and a mean of 36 calibration sources have been used in each frame. We have not found any 
appreciable color related trend. 

The calibrated images were inspected for the possible presence of ring patterns that 
would have been produced by pupil ghosts. We computed for each band, and for all the 
sources used in the calibration process of the different final frames, the difference between 
our photometry and the 2MASS photometry as a function of the radial distance to the 
nominal field center. 

We found a significant effect, over 0.02 mags only in the J-band images, as it is shown 
in Fig. HI This effect was c orrected by fi tting the pupil ghost, using the mscpupil task in the 



IRAF MSCRED package (jValdea 120021 ) . and removing this pupil image from the flat-fields 
and individual images. The final resulting radial differences are shown for the J band in 
Fig. [5] where it clearly appears that all the systematics was corrected well below the la level. 



3. Galaxy Number Counts 



The steps to co mpute robust counts to the 50% d etection level of the images are sim ilar 
to those described in ICristobal-Hornillos et al.l (120031 ) and in lEliche-Moral et al.l (120061 ) . In 
the next section we detail how the best set of SExtractor parameters was estimated as a 
compromise between optimizing the depth at the 50% completeness level while keeping low 
the number of spurious sources. 



3.1. Completeness corrections 

In order to compute the corrections that should be applied to the faint part of the 
galaxy counts we have performed a set of Monte Carlo simulations where real sources from the 
images were injected back to the same science image at different positions. The completeness 
correction to be applied depends on the surface brightness profile of the source. To account 
for this we computed a different correction function for sources in three effective radius (Re) 
intervals. Those Re intervals in pixels for the simulations were chosen from the histogram 
of the Fig. 1 (top panel): Re < 1.75, 1.75 < Re < 2.25 and Re > 2.25. In Fig. 1 (bottom 
panel) it is shown how the Re recovered decreases as the magnitude goes fainter (see below 
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Fig. 2. — Calibration with 2MASS stars for the one of the pointings included in the ALH08 
field in the H filter. Sources marked with a circle were eliminated using a 2a reject in the 
calibration fit. Error bars correspond to the 2MASS magnitude error summed in quadrature 
with the SExtractor computed error. 
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Fig. 3. — Histogram of the photometric accuracy of the different pointings in the ALH08 
field. In the small panel the histogram of the number of used stars for the calibration in 
each pointing is shown. 
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Fig. 4. — Differences between 2MASS and ALHAMBRA photometry versus the distance to 
the nominal pointing center in the J band. 16 pointings have been combined and a binning 
have been performed in the x-axis to increase the signal. For each bin the boxplot shows 
the mean (stars), the median, 1st and 3rd quartile (boxes), the minimum and maximum 
(circles), values after trimming the 10% on either side of the sample in each bin. The error 
bars represent the la interval. 
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Fig. 5. — Differences between 2MASS and ALHAMBRA photometry versus the distance 
to the nominal pointing center after applying the pupil ghost correction in the flat-field as 
explained in the text. Symbols as in Fig. HI Note that the scale in the y-axis has change 
with respect to Fig. HI 
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for a description of how these values were obtained). 

For the simulations, bright sources were selected from the image in each Re interval. 
These sources were artificially dimmed to the 0.5 magnitude bin under study and injected 
back on the image randomly. In each iteration 40 sources were simulated, computing the 
recovering fraction and the robust mean (trimming the 10% on either side of the distribution) 
of the SExtractor desired output parameters differences between the dimmed input sources 
and the recovered ones. Using 9 of these iterations in each magnitude bin, we estimated the 
mean and rms of the recovered fraction, and SExtractor parameters differences over all the 
meaningful magnitude range. In Fig. [TJ the completeness correction for the three Re ranges 
is shown for one of the pointings in the J band. 

It can be argued that with that method unrealistic pseudo-artificial sources could be 
produced. Whereas this could be true, the goal of this procedure is to parameterize the 
recovering efficiency on the basis of the source and image characteristics without any physical 
consideration, which will be implicitly taken into account when performing the corrections 
on the real data. As a validation of the procedure, the same simulations have been performed 
using real NICMOS F110W sources from the HDF-S Flanking Fields. The NICMOS images 
were resampled and convolved with a suitable Gaussian kernel to match the pixel scale and 
FWHM of the OMEGA2000 image under study. Initial sources were taken in the interval [m- 
l,m+l] (being m the magnitude under study), which produces a realistic mag — R e relation. 
As can be seen from the Fig. [7J the results in the completeness correction are the same that 
using dimmed sources from the image. 

In the top panel of Fig. [8]we show the depth to the 50% and 80% of recovering efficiency, 
computed using a linear interpolation in the magnitude vs. efficiency data. The figure 
indicates that a decreasing of the detection thresholds, in order to get a fainter level at 
the 50% of recovery efficiency limit, will not improve by much the limit at the 80% of 
completeness, which seems to reach a plateau. Moreover, as will be seen from the reliability 
plots, it produces a significant increase of detected spurious sources, as we explain below. 

3.2. Detection reliability 

To accurately compute the galaxy number counts it is important to establish the re- 
liability of the detections at the faint magnitude end. To find out the optimum way to 
evaluate that reliability we have studied the performance of three different methods. The 
first approach was to create artificial sky images with the same rms and background distri- 
bution than the real ones. The ratio of spurious versus real detections was computed running 
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Fig. 6. — (Top panel) Histogram of Re for the sources in one of the pointings of the field 
ALH08 (Bottom panel) Differences between Re(in) and Re(out) for the dimmed injected 
sources. The symbols correspond to x Re(in)< 1.75 pixels, □ 1.75 <Re(in)< 2.25 and A 
Re(in)> 2.25 pixels. The shaded is the range of magnitudes between the 50% and the 80% 
of detection efficiency. 
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Fig. 7. — Completeness correction for the same pointing of the previous figures in the J band 
using 0.8<7 DETECT_THRESH and the 2.5 pixel gaussian kernel. The vertical dotted (solid) 
lines are the magnitude at the 50% (80%) of detection efficiency. 
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Fig. 8. — (Top panel) 50% vs. 80% of recovering efficiency for the same pointing as before in 
Ks band using different SExtractor detection strategies. The size in pixels of the SExtractor 
Gaussian convolution kernel is given in the labels. The detection thresholds are indicated at 
the points. (Bottom panel) 50% of detection efficiency vs. spurious to total ratio. 
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SExtractor over the science and the artificial sky images. 



We have also inspected the performance of the method used in lCristobal-Hornillos et al. 



( 120031 ). Basically half exposure time images were created from two complementary sets of 
the data. The detections were performed in the total time image and the source fluxes 
measured in the half time images using the SExtractor double image mode for the same 
automatic apertures. Those created sources showing a magnitude difference greater than 3cr 
were considered spurious. 

The last method and the one that, at the end, produces the best results, consisted in 
constructing sky images using a similar combination procedure that the used to create the 
science images: combining the unregistered processed images with subtracted background 
using an artificial dither pattern. The major difficulty here was to remove the extended 
sources that could bias the sky even when doing trimmed mean (discarding 20% of the 
pixels at each side). We have confirmed that SExtractor does locate these smooth deviations 
over the sky rms when filtering is used. To avoid this we multiplied those sky images by -1 
and used them as real sky images. 

In the Fig. [9] we show a comparison between the spurious rate at the 50% of detection 
efficiency computed over the 'real sky' and the artificial sky. A good agreement between 
both methods is observed, indicating that the use of artificial images to compute the ratio of 
spurious to total detections is adequate. Being artificial images faster to construct we made 
use of them to estimate the detection reliability in each pointing. 

The bottom panel of Fig. [8] shows the magnitude reached at the 50% of detection ef- 
ficiency versus the spurious to total ratio at the same magnitude bin. From the figure 
we can see that in order to reach a deep 50% detection limiting magnitude, maintaining 
at the same time the number of spurious detections below 20%, the optimum SExtractor 
DETECT JTHRESH-FILTER combinations are: DETECT_THRESH=1.2 or 1.4 without fil- 
tering, or DETECT_THRESH=0.8 using a filtering with a Gaussian kernel of size similar 
to the image FWHM. In the top panel of Fig. [8] it is shown that those combinations reach 
roughly the same magnitude at the 80% of recovery efficiency. We have decided to use the 
latter filter-DETECT_THRESH combination because, as can be outlined from the bottom 
panel in Fig. [10J the differences between the input and the recovered AUTO magnitudes in 
the simulations are close to zero in all the magnitude range up to the magnitude at the 80% 
of recovery efficiency. Close to the magnitude at the 50% of completeness the recovered mag- 
nitude appears to be ~ 0.1 — 0.2 mag brighter than in the previous bin. This effect could 
be due to the fact than only those sources that suffer from noise brightening were found 
by SExtractor and used to compute the input-output magnitude differences, producing a 
bias towards brighter recovered objects. Following those results, the number counts in the 
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Fig. 9. — Ratio of spurious to total detections when artificial and real sky images (see text) 
were used. Different SExtractor DETECT_THRESH (er) and FILTER (F) combinations are 
shown. The grey area corresponds to a difference less than O.f. 
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following sections will be computed and corrected up to the magnitude of 80% of recovery 
efficiency for point-like sources, avoiding any possible systematics due to a magnitude shift 
at the 50% completeness bin. 

As pointed before, the photometry of the sources was obtained using MAG_AUTO. 
Simulations indicate no significant differences between simulated and recovered magnitudes 
(at the 80% recovery efficiency) for the used values of the SExtractor parameters FILTER 
and DETECT_THRESH. Using the same kind of simulations we computed the rms of the 
recovered MAG_AUTO values in each bin to characterize the photometric error, the results 
are shown in Fig. [TTJ 



3.3. Star-galaxy separation 



In ICristobal-Hornillos et al.1 (120031 ) it was shown that fainter than Ks=17.0 the correc- 
tion due to stars is < 0.06 dex. However, given the lower Galactic latitude of the ALH08 
field, a higher number of contaminating stars is expected. A correct star subtraction is 
relevant in the intermediate magnitude range. At bright magnitudes stars can be easily sep- 
arated from galaxies using a compactness criteria. In contrast, at intermediate magnitudes 
the star/galaxy separation is more demanding because many galaxies are barely resolved 
with small apparent sizes compared with the FWHM and pixel resolution, as can be seen 
from the diagram shown in Fig. [121 

The viability of star/galaxy separation using the SExtractor neural network has also 
been analyzed. For this purpose, using the same Monte Carlo method explained before, 
bright stars and galaxies from the images have been artificially dimmed to each magnitude 
bin and their SExtractor stellarity parameter recovered. In Fig. [13] it clearly appears how, in 
one final frame with FWHM=1.1, the input and recovered CLASS_STAR differences are less 
than 0.1 up to Ks=18.0. But, in the next bin Ks> 18.5, the CLASS_STAR of the dimmed 
stars and galaxies could lead to some misclassification. Nevertheless, as can be seen from 
the histogram in the bottom panel of Fig. [TBJ a non negligible number of objects start to 
populate the range from 0.4 to 0.8 in CLASS_STAR at magnitudes fainter than 17 in Ks, 
and the selection of the CLASS_STAR cut off might bias the star counts estimates. 

An al ternative way to per form the star/galaxy separation makes use simply of color-color 



diagrams. iHuang et al.l (119971 ) have established a reliable star/galaxy separation using the B- 
I vs I-K colors. Here we make use of the SDSS DR5 data for the ALH08 field to proceed with 
this separation using the g-r vs J-Ks colors shown in Fig. [TH The star counts are corrected 
by the ratio of the Sloan/Alhambra completeness factors in the corresponding filter shown in 
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MAG_AUTO(input) 

Fig. 10. — Differences between the input and the SExtractor recovered AUTO magnitudes 
in one of the pointings of ALH08 in the J band for two SExtractor FILTER and DE- 
TECT_THRESH combinations: no filtering and DETECT_THRESH=1.2 (upper panel), and 
2.5 pixel filtering and DETECT_THRESH=0.8 (bottom panel). The vertical dotted (solid) 
lines are the magnitude at the 50% (80%) of detection efficiency. 
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Fig. 11. — Median photometric random errors in magnitudes per magnitude bin (<r m ) in 
the three NIR filters for the eight ALH08 pointings (crosses and dotted lines), the three 
lines represent sources with Re <= 1.75, 1.75 < Re <= 2.25, and Re > 2.25 pixels. The 
magnitude errors are computed as the rms of the recovered sources magnitudes in each bin 
of the simulations described in the text. Error bars represent the rms of the computed a m 
among the pointings. Exponential grow fit (a(m) = <j + aexp [b(m — m )]) to the magnitude 
errors in each band (solid lines). 
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Fig. 12. — Peak/ISOPHOTAL_AREA vs. Ks diagram. The red dots are object for which 
the SExtractor stellarity index is >0.8. 
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Fig. 13. — (Top panel) CLASS_STAR input-output differences as a function of input artifi- 
cially dimmed magnitude for a set of well defined stars (CLASS_STAR(in) ~ 1) and galaxies 
(CLASS_STAR(in) ~ 0). In the case of stars l-(CLASS_STAR(out)-CLASS_STAR(in)) is 
plotted. (Bottom panel) Histogram of CLASS_STAR for different ranges in the Ks magni- 
tudes. 
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Fig-HSl The star counts using Sloan-Alhambra colors were computed to magnitude 19.5,18.5 
and 18.0 respectively in the J,H and Ks band, where the Sloan/ Alhambra completeness is 
>0.5. For the fainter star counts we used a 0-slope extrapolation where the models of star 
counts in the Galaxy are flatter (see Fig. [TBI) . In Fig. [15] bottom panel, the correction to the 
log(iV) galaxy counts is presented, showing that fainter than the magnitudes up to where 
the color-color separation could be performed the correction in the three bands is < 0.08 
dex. 

To check the validity of the sta r counts computed using the color-color approach, in 



Fig. US] those are compared with the iRobin et al.l (119961 ) Galactic star counts models. The 



lines correspond to star counts in the Galaxy for galactic coordinates (l=100,b=-45), very 
close to the coordinates of the Alhambra-08 field (1=99, b=-44). From this figure it appears 
that it is possible to accurately remove the stars from the galaxy counts using the color-color 
diagram. This was the method we have finally used to do the star/galaxy separation. 

Finally is noteworthy to mention that, star/galaxy separation will be done more ac- 
curately when the photometry in the 23 Alhambra bands is available, as the system will 
provide a sort of low resolution spectroscopy for each object. 



4. Results 

The corrected galaxy number counts have been computed for the ALH08 field in the 
three standard NIR filters in a consistent way in the sense that they have been estimated 
following the same scheme for the three bands. The number count data, computed and cor- 
rected up to the magnitude of 80% of recovery efficiency for point-like sources, are presented 
in the Tabs. H d and [3j 

The error in the number counts for each pointings are the sum in quadrature of the rms 
in the estimation of the completeness corrections, with the contribut ion of Poisson noi se and 
galaxy clustering calculated for each magnitude bin following eq. [Tl (jliuang et al. 1997 ). that 
includes the angular correlation function. 

a\ = Ni(m) + 5.3 ^ )/2 N^m) (1) 



being iVj the raw counts in the pointing, r = 7.1 Mpc, 7 = 1.77 and 51og(r*) = 
m - M* - 25. M* is set to -22.95,-23.69 and -23.93 for the J, H and Ks filters. The final 
uncertainty in the combined counts per square degree and magnitude is given by: 
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J-Ks(AB) 



Fig. 14. — g-r vs J-Ks plot used to perform the star-galaxy separation. All objects in ALH08 
with S/N > 5 are plotted. The line separating the two classes is g — r = 1.2 + 2.0-(J—Ks)ab- 
The plotted symbols correspond to the Sloan classification dots for representing stars and 
the crosses for galaxies. R edshift-color tracks for three galaxy models constructed using the 



Bruzual fc Charlotl (120031 ) code are shown. The redshift values are displayed over the Single 



Stellar Population (SSP) track, the marks over the other 2 cur yes have the s ame r edshift 
spacing. The diamonds represent the mid stellar locus given in ICovey et al.l (120071 ). Two 
contour lines containing the 95% and the 68% are displayed for the galaxies and stars using 
the Sloan classification. 
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Fig. 15. — (Top panel) The ratio of sources identified in the Alhambra field for which there 
is a counterpart in the Sloan SDSS DR5 catalog. (Bottom panel) The correction subtracted 
from log(N) galaxy counts due to stars. The big open circle indicates the point up to which 
the number of stars can be estimated from the Sloan- Alhambra colors (see text), the fainter 
star counts assume a flat extrapolation from this value. 
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Fig. 16. — Star counts for the ALH08 field compared with the iRobin et al.l (119961 ) Galactic 
count model. The J and Ks data have been plotted with an offset of -0.5 and 0.5 dex in the 
y-axis. 
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In the Figs. [13 [TBI an d dH the Alhambra counts in the J,H, and Ks bands are plotted 
together with the computed number counts from other surveys. 

Let us first point out the general aspect of those counts, leaving more detailed consider- 
ations for the next section. Our c orrected J band galaxy counts are in ove rall agreement with 
those computed in othe r works (ITeplitz et al.lll999l ; iFeulner et al.l 120071 ). t he bright end of 



the re sults presented by lMaihara et al.l (120011 ) . and the ELAIS-N2 data from I Vaisanen et al. 



( 120001 ) . Nevertheless, our counts are lower than those given by iBershady et al.l (119981 ) at 
magnitudes J > 21.5. 

The Alhambr a counts in the H ba nd are in good correspondence with the pub lished data 



by iMartinil (120011 ); iFrith et al.l (120061 ) and with the bright part of the data from iMoy et al. 



d2003h . after applying an offset of -0.215 mag. This offset was calculated following th e -0.06 5 
calibration difference among Las Campanas Infrared Survey (LCIRS, Chen et al. ( 20021 )) 
reported in IMoy et al.l (120031 ). and the syste matic of -0. 2 8 ma gnitude difference between 
LCIRS and 2MASS magnitudes reported in IFrith et al.l (120061 ). However, the faint end 



(H > 20) of our data is s ignific antly above the faint number count data from lYan et al. 
fjl998h and llMetcalfe et al.l (bood ). obtained from NICMOS observations. 



Regarding the Ks filter, for which there are numerous number counts studies, our galaxy 
counts are in good agreement with most of the published data, as can be appreciated in the 
Fig. | 



4.1. The measured slopes 

We have found, as in previous studies, that the slope of the galaxy number counts 
displays a clear change at Ks~17.3 (see Fig. [2D]). However, in the present work we also have 
found this change of slope in the J and H band galaxy counts at J~19.1 and H~ 18.2. 

The slopes of the bright end and faint end counts in the different bands were measured 
using the least squares method. The values that we have calculated and the magnitude 
ranges used for the fits are listed in Tab. [H The results show that the faint slopes for the 
three NIR filters are the same within 2a, whereas the bright slope is steeper in the Ks band. 



As is described in iTeerikorpil ( 120041 ) the fact that the photometric error increases at 
fainter magnitudes, and that the differential galaxy number counts also rises towards lower 
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Table 1. Corrected galaxy number counts in the J band. 



Magnitude raw counts a eff. cor. b log(iVc) log(Nm) d area 

N ■ mag~ l • deg~ 2 N • mag^ 1 ■ deg~ 2 deg 2 



13.00 


43.0 


1.00 


1.31+?-*? 

— l.ol 


1.30+°^ 


0.441 


13.50 
14.00 
14.50 


31.0 
53.0 
75.0 


1.00 
1.00 
1.00 


92+ - 73 
u - yz -0.92 

1 iq+0.66 

^^-l.lS 
1 63+ - 36 


o.95l8:^ 
Losing 


0.441 
0.441 
0.441 


15.00 
15.50 


118.0 
129.0 


1.00 
1.00 


i ro+0.28 
1.88_ 107 

z - uu -0.54 


1 aq+0.49 
1 - oa -1.89 

2 00+ - 41 
z.uu_ 2 00 


0.441 
0.441 


16.00 


175.0 


1.00 


o 11+0.22 
z - ii -0.47 


9 i i +0.34 
z - ii -2.11 


0.441 


16.50 


242.0 


1.00 


o ec+0.10 
z - oo -0.13 


Z - OD -0.55 


0.441 


17.00 


309.0 


1.00 


9 77 +0.07 
z -' '-0.08 


9 77+O.I5 
z -' '-0.23 


0.441 


17.50 
18.00 


436.0 
634.0 


1.00 
1.00 


9 qi+0.06 
z - yl -0.07 
o 99+O.O4 
°' zz -0.04 


9 Qi +0.09 
z - yl -0.12 
q 99+0.11 
°- zz -0.14 


0.441 
0.441 


18.50 


832.0 


1.00 


o 4Q+0.03 
°- 4U -0.03 


o 4f) +0.09 
°- 4U -0.11 


0.441 


19.00 


1128.0 


1.07 


O R q + 0.02 

°- uo -0.02 


o 64 +0.06 
°- D4 -0.07 


0.441 


19.50 


1615.0 


1.08 


q on+0-02 
^• 8U -0.02 


q on+0.05 
d - 8U -0.05 


0.441 


20.00 
20.50 
21.00 


2303.0 
3285.0 
4455.0 


1.09 
1.12 
1.17 


Q, qq+0.06 
°- yy -0.08 

4-1818:81 
4.34l8;8i 


Q, qq+0.08 

°- yy -o.io 

A iq+0.06 
4 - i8 -0.07 
4 Q4+0.05 
4 - 34 -0.05 


0.441 
0.441 
0.441 


21.50 
22.00 


5030.0 
1739.0 


1.25 
1.52 


4 - ou -0.03 

4 «7+o.03 

4 - ' -0.03 


4 ^n+ 03 

4 - ou -0.03 

4 67+ - 01 

4 - D '-0.02 


0.381 
0.109 



a Raw counts including stars. 

b Effective correction, defined as (Nc + Nst)/(2 ■ area ■ raw). 

c Corrected galaxy counts,errors corresponds to the Poissonian and galaxy clustering 
uncertainty plus error in completeness added in quadrature. 

d Mean of the 8 pointings and its rms. 
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Table 2. Corrected galaxy number counts in the H band. 



Magnitude 


raw counts a 


eff. cor. 


log (iVc) 

AT —1 7 —2 

TV • mag ■ deg 


log(iVm) a 
iv • mag ■ deg 


circci 
deg 1 


12.00 


22.0 


1.00 


97 +0 n 6 Q i 

' —0.97 


98+2-2! 


0.444 


12.50 


40.0 


1.00 


1 62+°i° ) 

u — 1.62 


1 60+?-^ 

- L-uu — 1.60 


0.444 


1 3 00 


39 


1 00 


1 iq+0-56 
i - iy -1.19 


1 TO+0.79 
llly -1.19 


0.444 


1 3 50 


53 


1 00 

X . \J\J 


52+ 119 


55 +L44 


0.444 


14 no 


Q7 


1 DO 


1 85+ - 28 


1 qr+0.49 
A - 80 -1.86 


n 444 


14.50 


106.0 


1.00 


1 04+0.24 


1 94+0.54 
± - y ^-1.94 


0.444 


15.00 


131.0 


1.00 


1.8ll?-jj? 


i.8ii?i? 


0.444 


15.50 


223.0 


1.00 


z.oy_ 21 


9 qn+0.25 
Z - Oy -0.63 


0.444 


16.00 


288.0 


1.00 


2 74 +o.07 

z - '^-0.09 


9 74 +o.i6 

z - '^-0.25 


0.444 


16.50 


363.0 


1.00 


O OO+0.06 


9 OO+0.12 


0.444 


17.00 


541.0 


1.00 


o in+0.04 
°- iu -0.05 


ln +0.06 


0.444 


17.50 


788.0 


1.00 


o O7+0.03 
-0.03 


q q 7 + 0.10 

°-°'-0.14 


0.444 


18.00 


1074.0 


1.00 


o r; fi +0.02 
J - oo -0.02 


q t;7+0.06 
o - d '-0.07 


0.444 


18.50 


1512.0 


1.05 


o 7fi +0.02 
°- ' u -0.02 


q 7fi +0.06 
°- ' D -0.07 


0.444 


19.00 


2110.0 


1.05 


o Q4+0.07 
°- y ^-0.08 


04+O.O8 


0.444 


19.50 


3138.0 


1.07 


4 14+0.04 

4 ' 14 -0.05 


4 14 +o.05 

4 - i4 -0.06 


0.444 


20.00 


4456.0 


1.12 


A qn+0.03 
^• oz -0.03 




0.444 


20.50 


5879.0 


1.23 


4 40+0.02 

^•^-0.02 


4 40+0.04 

4 - 4y -0.05 


0.444 


21.00 


2578.0 


1.48 


4 - DO -0.02 


4 -°°-0.02 


0.165 



a Raw counts including stars. 

b Effective correction, defined as (Nc + Nst)/(2 ■ area ■ raw). 

c Corrected galaxy counts,errors corresponds to the Poissonian and galaxy clustering 
uncertainty plus error in completeness added in quadrature. 

d Mean of the 8 pointings and its rms. 
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Table 3. Corrected galaxy number counts in the Ks band. 



Magnitude 


raw counts 


rr V) 

en. cor. 


log (TV c) 

AT —1 J —2 

TV • mag ■ cleg 


1 / A7 \ d 

log(iVmJ 

AT —1 J —2 


area 
deg z 


12.00 


25.0 


1.00 


J --- L -1.14 


1 13+? 


0.440 


19 50 


44 


1 00 


1 61+ - 32 


1 60+ - 52 
1 - ou -1.60 


440 


lo.UU 


a i n 

41. U 


1 nn 

l.UU 






n 4 /in 

U.44U 


13.50 


78.0 


1.00 


i"-89±8:it 




0.440 


14.00 


101.0 


1.00 


9 no+0.20 
Z - UO -0.39 




0.440 


14.50 


141.0 


1.00 


9 91+O.I6 
z - zl -0.26 


9 99+O.32 
Z.z:z_2.22 


0.440 


15.00 


188.0 


1.00 


2 41+0.12 
z - 4i -0.17 


2 41 +0.26 
z - 41 -0.73 


0.440 


15.50 


291.0 


1.00 


2 74+0.07 

z - ' 4 -0.09 


2 74+°- 14 

z - ' 4 -0.21 


0.440 


16.00 


375.0 


1.00 


z - yu -0.06 


2 96 +an 
z - yD -0.15 


0.440 


16.50 


585.0 


1.00 


o 99+O.O3 
°' zz -0.04 


q 99+O.O8 
°- zz -0.09 


0.440 


17.00 


940.0 


1.00 


o 40+O.O2 
°- 4y -0.02 


40+O.O6 
°- 4y -0.07 


0.440 


17.50 


1307.0 


1.03 


O C7+0.02 

°- u '-0.02 


q C7+0.05 
°- u '-0.06 


0.440 


18.00 


1836.0 


1.04 


3.86ig:S 


3-8618:8? 


0.440 


18.50 


2614.0 


1.05 


4-o4±g:8i 




0.440 


19.00 


3698.0 


1.11 


4 94+0.04 
4 - Z4 -0.04 


4 24+ - 06 
4 - Z4 -0.07 


0.440 


19.50 


4477.0 


1.37 


4 49+0.03 

4 - 4Z -0.03 


4 - 4z -0.04 


0.440 


20.00 


574.0 


1.55 


4 ^n+ - 04 

4 - ou -0.05 


4 ^n+ - 02 

4 - ou -0.02 


0.054 



a Raw counts including stars. 

b Effective correction, defined as (Nc + Nst)/(2 ■ area ■ raw). 

c Corrected galaxy counts,errors corresponds to the Poissonian and galaxy clustering 
uncertainty plus error in completeness added in quadrature. 

d Mean of the 8 pointings and its rms. 
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J(Vega) 



Fig. 17. — Galaxy number counts in the J filter com pared with data from other sur veys. 
The lines correspond to t he number counts models in ICristobal-Hornillos et al.l (120031 ) and 
Eliche-Moral et al.l (120061 ) described in the text. 
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Fig. 18. — Galax y number counts in the H filter com pared with data from oth er surveys 
and the models in ICristobal-Hornillos et al.l (120031 ) and Eliche-Moral et al.l (120061 ) described 
in the text. 
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Fig. 19. — Galaxy number coun ts in the Ks filter compared with data from the literature 
and the galaxy counts models in ICristobal-Hornillos et al.l (120031 ) lEliche-Moral et al.l (120061 ) 
described in the text. Also it is shown a model where only the passive evolution of stellar 
populations is considered (PLE). 
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Intercept: -2.852 ±0.203, Slope: 0.358 ±0.010 
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Fig. 20. — Galaxy number count bright and faint slopes found in this work. Error bars are 
the Poissonian and galaxy clustering uncertainty added in quadrature with the rms in the 
estimation of the completeness corrections. 



-36- 



fluxes, lead to a stepper observe d slope at the faint end. This effect is related with the 
Eddington bias (lEddingtonl Il913l ) . To investigate if the computed magnitude errors could 
bias our slope estimates we have done the following study. First, in each filter we use 
the exponential grow function fit to the magnitude errors (Fig. [TTT) . Then we extended 
our computed number counts one magnitude fainter using the corresponding slope in Tab. HJ 
simulating a Gaussian decay for fainter magnitudes. The parameter a for the Gaussian is the 
value of the exponential fit to the error values for point-like sources at one magnitude fainter 
that the last bin given in Tabs. [TJ [2] and [3] (er = 0.35,0.41 and 0.44 at respectively J=23.0, 
H=22.0 and Ks=21.0). Using the fitted magnitude errors for point-like objects (cr(ra)), and 
assuming that real distribution is close to the extended number counts (N(m)), we simulate 
the bias produced by the photometric error on the observed counts, by convolving the N(m) 
with cr(m) using the eq. The results indicate an increase in the slope less than 0.015 in 
the case of J and H filters, and 0.03 in the Ks band, meaning that the original distributions 
would have a faint end slope of 0.33,0.34 and 0.30 in the J,H, and Ks filters in the ranges 
given in Tab. HI 



N(m. 



obs I 



N(m) 



V2 



(m-m ob3 ) 

-e 2ct2 (™) dm 



no 



(3) 



The increase of the slope is not observed when we simulate this bias directly over the 
observed count distribution. This is due to the fact that fainter than the 80% completeness 
magnitude bin there is a fast decrease of the number of detected sources, and that in this 
bin the typical photometric error [o m = 0.23 — 0.26) for the dominant point-like sources, 
makes that the bias do not significantly affect to the previous bins. In this case the slopes 
after applying the bias are lower than the computed from the observed number counts 
(0.33 ±0.01, 0.33 ±0.02, and 0.30 ±0.02 in J,H and Ks). This result would suggest that the 
original distribution slopes at the faint end would be higher than the ones given in Tab. HI 
We have used the combined results from this paragraph and the previous one to increase 
the uncertainty in the slope, leaving the values computed directly from the observed number 
counts as satisfactory estimates of the real distribution slope at the faint end. 

Our results for the Ks band are in good agreement with other K-band surveys which also 



report a simi l ar change of slope in the galaxy counts in the range 17.0-1 8.0 (IDaddi et al.ll2000 



Huang et al.l l200ll ; ICristobal-Hornillos et al.l l2003l; llovino et al.l I2005T). A t the bright par t 



our slopes are close to the valu es measured in 



Cristobal-Hornillos et all (J2003J) 



point to a steeper bright slope (I Gardner et al 



Iovino et al 



Kiimmel fc Wagner! fl2000[ ); Martini (120011 1; 
(2005h see Tab. ED, while other ref erences 



19931 . Il996l : Huang et al.lll997L 1200 if ). This 



could be due to the fact that the last authors could extend the power-law fit to brighter 
magnitudes due to their larger surveyed area, whereas our fit is closer to the magnitude 



-37- 



where the break is found which would lead to a decrease in the slope if the break transition 
is smooth. In the faint part o f the K s counts a s l ope of 0.33 found in thi s work is in 
agreement with iBershady et al.l (119981 ); iHuang et al.l (120011 ); llmai et al.l (120071 ). The value 
of the faint slope from the ALHAMBRA data is however steeper than the value reported in 
some surveys coverin g smaller areas, where th e fitted range ex t ends to fainter magnitudes, 
Gardner et al.l (119931 ); iMoustakas et al.l (119971 ); iMaihara et al.l ( 120011 ). Here the differences 
might b e due to cosmic variance . A la rger than 2a disagreement is found with the faint 
slope of ICristobal-Hornillos et al.l (120031 ) fitted in the range [17. 5,19.51, although thei r value 
increase to 0.29 when the fit interval is extended to Ks=21.0. iKummel fc Wagnerl (120001 ) 
give a higher value for the faint slope, however due to the brighter limiting magnitude of 
their number counts they could established a break or the beginning of a roll-over in the 
interval Ks=[16.5,17.0]. 

In the H and J bands there are fewer works reporting count-slope values. Our results, 
given in Tab. HI show similar slopes for the J and H filter at the bright and faint ends. As 
can be seen in Tabs. [6] and [7] our result at the b r ight e nd a re in good agreeme nt with the 
bright-end slope values to H=19 given in 
J filter with the results in 
only the slope va 



ues m 



Vaisanen et al. 



Martinil (120011) andlChen et al l (120021 ) . and in the 



jgoooh 



and 



Iovino et al.l (120051 ) . At the faint end, 
Maihara et al] ( 2001 ) in the J filter, estimated in an area of 4 arcsec 2 , 



and I Chen et al.l (120021 ) have a significant discrepancy. 



5. Comparison with Models of Evolution 

Historically, the galaxy number counts have been used to examine parameters of the 
cosmological model and to test different galaxy evolution scenarios. Now that the cosmo- 
logical parameters are fixed using other methods the consequences derived from the galaxy 
counts for the galaxy evolution have become more precise. The ALHAMBRA computed 
counts in the three NIR bands provide a good dataset which, when combined with other op- 



Table 4. Measured slopes in the J,H and Ks filters 



Filter 


Bright range 


Bright slope 


Faint range 


Faint slope 


J 


[17.0,18.5] 


0.44±0.04 


[19.5,22.0] 


0.34±0.01 


H 


[15.5,18.0] 


0.46±0.02 


[19.0,21.0] 


0.36±0.02 


Ks 


[15.0,17.0] 


0.53±0.02 


[18.0,20.0] 


0.33±0.03 
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Table 5. Characteristics of the surveys in the K band 



Reference 


Surveyed area 
sqarcmin 


limit magnitude 


bright range 


bright slope 


faint range 


faint slope 


filter 


Gardncr93 


5688 


14. 5 C 


[10.0,16.0] 


0.67 


[18.0,23.0] 


0.23 


K' 


Garrinor93 


582 


16.75 c 












VJJCX1 U11C1 kJ*-J 


167.7 


18.75 c 












G&rdricr93 


16.5 


22.5 












Glazcbrook94 


552 


16. 5 C 










K 


1~) inrPTiv^Wi Q.'S 


3 


23. 5 C 






[20 23 51 


32+0 02 


K 


MrLeodQR 


22.5,2.0 


19 5 21 25 c 










Ks 


Gardner96 


35424 


15.75 c 


<16.0 


0.63±0.01 


— 


— 


K 


Moustakas97 


2.0 


24.0 C 






[18.0,23.0] 


0.23±0.02 


K 


Huang97 


29628 


16. C 


[12.0,16.0] 


0.689±0.013 






K' 


Minezaki98 


181,2.21 


19.1,21. 2 d 






[18.25,18-75] 


0.28 


K' 


Bcrshady98 


1.5 


24.00 b 






>18.5 


0.36 


K 


Szokoly98 


2185 


16.5 C 


[14.5,16.5] 


0.50±0.03 






Ks 


Saracco99 


20 


22.25 c 






[17.25,22.5] 


0.38 


Ks 


VaisancnOO 


3492,2088 


16.75,17.75 c 


[15.0,18.0] 


0.40-0.45 






K 


MartiniOl 


180,51 


17.0,18.0 C 


[14.0,18.0] 


0.54 






K 


DaddiOO 


701,447 


18.5,19.0 C 


[14.0,17.5] 


0.53±0.02 


>17.5 


0.32±0.02 


Ks 


KiimmclOO 


3348 


17.25 c 


[10.5,17.0] 


0.56±0.01 


[16.5,17.5] 


0.41 


K 


MaiharaOl 


4 


25.25 c 






>20.1 


0.23 


K' 


SaraccoOl 


13.6 


22.75 c 






>19 


0.28 


Ks 


HuangOl 


720 


19. 5 C 


<16.5 


0.64 


>17 


0.36 


K' 


Cimatti02 


52 


19.75 c 










Ks 


Cristobal03 


180,50 


20.0,21. b 


[15.5,17.5] 


0.54 


[17.5,19.5] 


0.25 


Ks 


Iovino05 


414 


20.75 c 


[15.75,18.0] 


0.47±0.23 


[18.0,21.25] 


0.29±0.08 


Ks 


Elston06 


25560 


19. 2 b 










Ks 


Imai07 


750,306 


18.625,19.375 c 


<18.00 


0.32±0.06 


>19.00 


0.32±0.06 


Ks 


Feulner07 


925 


20.75 c 










K 


This work 


1584,194 


19.5,20.0 C 


[15.0,17.0] 


0.53±0.02 


[18.0,20.0] 


0.33±0.03 


Ks 



a 3<r limit for point sources 

b 50% efficiency for point objects 

c the latest magnitude bin in the number counts 

d 80% efficiency for point objects 
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Table 6. Characteristics of the surveys in the J band 



Reference 


Surveyed area 
sqarcmin 


limit magnitude 


bright range 


bright slope 


faint range 


faint slope 


filter 


Bershady98 


1.5 


24.5 b 






>19.5 


0.35 


J 


Teplitz99 


180 


21.75 c 










J 


Saracco99 


20 


23.75 c 






[18.0,24.0] 


0.36 


J 


VaisancnOO 


2520,1275 


18.25,19.25 c 


[17.0,19.5] 


0.40-0.45 






J 


MartiniOl 


180,27 


18.5,20.5 C 


[16.0,20.5] 


0.54 






J 


MaiharaOl 


4 


26.25 c 






[21.1,25.1] 


0.23 


J 


SaraccoOl 


13.6 


24.25 c 






>20 


0.34 


J 


Iovino05 


391 


22.25 c 


[17.25,22.25] 


0.39±0.06 






J 


Feulner07 


925 


22.25 c 










J 


Imai07 


750,306 


19.625,20.375° 


[17.0,19.5] 


0.39±0.02 


>19.5 


0.30±0.03 


J 


This work 


1588,392 


21.0,22.0 C 


[17.0,18.5] 


0.44±0.04 


[19.5,22.0] 


0.34±0.01 


J 



a 3<r limit for point sources 

b 50% efficiency for point objects 

c the latest magnitude bin in the number counts 



Table 7. Characteristics of the surveys in the H band 



Reference 


Surveyed area 
sqarcmin 


limit magnitude 


bright range 


bright slope 


faint range 


faint slope 


filter 


Teplitz98 


35.4 


22. 8 a 










F160W 


Yan98 


8.7,2.9 


23.5,24.5 C 






[20,24.5] 


0.315±0.02 


F160W 


Thompson99 


0.7 


27.4 b 










F160W 


MartiniOl 


180,80 


18.0,19.0 C 


<19 


0.47 






H 


Chen02 


1408 


20. 8 d 


<19 


0.45±0.01 


>19 


0.27±0.01 


H 


Moy03 


97.2,619 


20.5,19.8 b 










H 


Metcalfe06 


49 


22. 9 a 










H 


Metcalfe06-Nicmos 


0.90 


27.2 C 










F160W 


Frith06 


1080 


17.75 c 










H 


This work 


1598,594 


20.5,21.0° 


[15.5,18.0] 


0.46±0.02 


[19.0,21.0] 


0.36±0.01 


H 



a 3cr limit for point sources 

b 50% efficiency for point objects 

c the latest magnitude bin in the number counts 

d S/N=5 in a 4" diameter aperture 
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tical data and independent determinations of the local luminosity functions, allow evolution 
to be examined, in particular the still uncertain question of the formation and evolutionary 
history of early type galaxies. 

In this section we compare our c ounts with semi -analytic predictions from number-count 
models, following the recipes given in lGardnerl (119981 ). which trace back the redshift evolution 
of the galaxy Spectral Energy Distri bution (SEP) of dif f erent g alaxy classes. The SEDs have 
been co mputed using the codes of iBruzual fc Charlotl (120031 ) . We apply dust attenuation 
following lEliche- Moral et al.l (120061 ) . r B = 0.6 that corresponds to Ty = 0.4 if the attenuation 
follows a oc A -2 power la w. We apply dust ex t inctio n either directly and by the same amount 
to all the galaxies using IBruzual fc Charlotl (120031 ) codes, following the prescription given 



in 



Chariot fc Fall! (120001 ). or by using the luminosity dependent extinction law proposed in 



Wangi (Il99ll ). The parameters we use to characterize four different galaxy types are given 
in Tab. In the LF parameterization for the different bands, M* is changed according to 
the rest-frame colors of the evolved SED (from Zf to z = 0), whereas a and 0* are assumed 
to be the same in all filters. 

In the first step we compare t he ALHAMBRA NIR counts w ith the prediction ob- 
tained using the model proposed in ICristobal-Hornillos et al.l (120031 ). The extinction cor- 
rection was applied direc tly to the SEDs, as an entry parameter in the code described in 
Bruzual fc Charlotl (120031 ) using Ty = 1.2 for stars younger than 10 7 yr , and fi = 0.3 as the 
fraction of it coming from an ambient contribution which affects the old stars too. The pa- 
rameters used in the local luminosit y function are M* = —24.07, a = —1.00, <ft* = 4.94 x 10~ 3 
calculated in iGardner et al.l (119971 ) (jCole et al.l (120011 ) provide the parameters for the A- 
cosmology), which were transformed to take into account the presence of different galactic 
types in the local LF adopting the galaxy mixing E/S0=28%, Sab/Sbc=47%, Scd=13%. 

The model also adds a dwarf star-forming population, characterized by an stellar pop- 
ulation of age lGyr at all redshifts, and a steeper slope LF (M* = —23.12, a = —1.5, 



0.96 x 10 3 ) given in 



Gardner! (119981 ). The formation redshifts are Zf = 2.0 for the 
E/S0 and intermediate-type disk galaxies and Zf = 1.0 for the Scd, although the forma- 
tion redshift for the Scd could be Zf = 4.0 or higher without modifying the total counts 
in NIR, that at magnitudes fainter than Ks=21.5 are dominated by the dwarf star- forming 
population. Due to the disappearance of the red-galaxy population at z > 2.0 this model 
reproduces the change of slope in the Ks band observed in the present data (see Fig. IT9|) . 



However, as it was discussed in lEliche- Moral et al.l (120061 ) . this model fails to simultane- 
ously reproduce the blue band counts as can be seen in Fig. EH where the predicted counts 
are compared with some B-band galaxy counts from the literature. 
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B(Vega) 



Fig. 21. — Galaxy number count in B band taken fro m the literature. The lines corres pond 
to the number cou n ts pre dictions from the models in ICristobal-Hornillos et al.l (120031 ) and 
Eliche-Moral et al.l (120061 ). Also it is shown a model where only the passive evolution of 



stellar populations is considered (PLE). 
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The number counts models that we consider now, below, include some modifications 
which aim to simultaneously reproduce the counts in both the NIR filters and blue filters. 
For these, the lo c al typ e dependent luminosity functions were computed from Sloan data in 



Nakamura et al.l (120031 ). as shown in Tab. O 



We consider first the two models proposed in lEliche-Moral et al.l (120061 ). In the first 
model a 0* evolution oc (1 + z) 2 , driven via mergers, is considered for the spiral and irreg- 
ular galaxies. This evolution in 0* is compensated by the evolution in M* to conserve the 
luminosity density. Is it important to take in mind that these models calculate the galaxy 
number counts tracing back the evolution of the stellar populations to z — 0, so the intrinsic 
brightening with z of the stellar populations must be added to the M* evolution when is 
compared with luminosity functions computed at higher redshifts. The formation redshift 
for the majority of the ellipticals and intermediate type disk galaxies in this model was set 
to 1.5. In the second model, the low formation redshift for the early spiral galaxies could be 
set to Zf = 4, avoiding at the same time an unreasonable high number of late type galax- 
ies at high-z, using the merger parameterization 0* oc exp [(— Q//3)((l + z)~^ — 1)] given in 



Broadhurst et al.l (119921 ). The value of = 1 + (2q n P 6 /2 was set to 1.53 using q = -0.55 



A value of Q = 1 was used as in lEliche- Moral et al.l (120061 ). The extinction correction, which 
i s im portant in the blue bands, was set to tb = 0.6 with the prescription given in IWang 



(119911 ). As is shown in Figs. [17] and HU these two models, that fit the B (Fig. [21]) and Ks 



(Fig. [19]) galaxy counts, overestimate the slope variation at J~19 in the Alhambra counts, 
and at H~18 in data from other surveys. 

In order to explain the change of slope in the NIR galaxy counts, the population of red 
Elliptical galaxies has to decrease with the redshift. Although a model taking into account 
only the stellar evolution with look back time fits the blue band counts (see Fig. [21]) . this 
model over-predicts the faint counts in the NIR bands as can be seen from Fig. [JJB Due 
to the red color of the slope change only the Elliptical population parameterized with a 
short burst of star-formation can play this role. Fig. [22] shows the evolution with redshift 

Table 8. Parameters for the SEDs 



Galaxy type functional form r Z/Zq IMF 

E/S0 Single star pop. 1 Salpeter 

early Sp Exponential 4 1 Salpeter 

late Sp Exponential 7 2/5 Salpeter 

Im Constant .. 1/5 Salpeter 
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of the J-H and J-Ks colors for an Elliptical galaxy and an early spiral formed at z=4 (with 
stellar populations accor ding to the parame t ers in Tab. El and reddening in the spectra 
applied directly from the iBruzual &: Charlotl (120031 ) code using ry = 0.8 and fi = 0.5). The 
estimated colors of the slope change from the fits given in Tab. [Hare J-H~ 0.97 ± 0.03 and 
J-Ks~ 1.84 ± 0.03, providing evidence that the Elliptical population at z~l must be the 
responsible of the turn down of the NIR counts. 

In the next model we introduce number- density evolution of the elliptical population 
parameterized using 0* oc (1 + z)~ 2 . The formation redshifts is set to Zf = 4.0 for all 
galaxy types, being in more agreement w ith the evolved red galaxies found at z > 2 
( Ivan Dokkum et al.ll2003t iDaddi et al.ll2004l ). The evolution in the early type densities was 
not accompanied by an evolution in M*, arguing that a substantial number of ellipticals 
formed in spiral-spiral mergers as expected for hierarchical galaxy formation. For the early 
spiral galaxies no number-density ev o lution was considered, and the density evolution pa- 
rameterization in lEliche-Moral et al.l (120061 ) 0* oc (1 + z) 2 was used for the two later type 



galaxi es. The reddening in the spectra was applied directly from the IBruzual fc Chariot 
( 120031 ) code using ry = 0.8 and fi = 0.5. This model, as can be seen in the Figs. [23], EH [251 
and EH] fits the galaxy counts in the optical and NIR filters, reproducing the feature of the 
slope turn down in the three NIR filters. 

In order to avoid an unreasonable high number of late ty pe galaxies at high - z, the simple 
merger evolution 0* oc exp [(— Q/ (3)((1 + z)~^ — 1)] given in iBroadhurst et al.l (119921 ) could 
be used with similar results. In this case, we used Q=-3 to parameterize the decrease in 
number of E/S0 galaxies, and Q=l to produce the required merger rate in the late spirals 
and irregulars. The number evolution given by these parameterizations are displayed in 
Fig [271 showing that at z ~ 1 the implied number density of E/S0 galaxies is only ~ 1/4 
of the present day 0*(O). Density evolution in the early type galaxies was observed in pre 



vious works studying the type-dependent LF evolution (iKauffmann et a 



2001 



Aeuerri fc Truiillol 12002 ; IWolf et al.l 120031 : iGiallongo et all 120051 : 



1996; 



Fried et al. 



Ilbert et al. 



2006 . 



Wolf et al.l (120031 ) found an increase in 0* of an order of magnitude for the early type galax- 
ies from z ~ 1.2 to z — 0, that is over the 0* oc (1 + z)~ 2 simulated here. However, they 
used the spectra of a present day Sa type galaxy to separate the different galaxies, which 
leads to an over-estimation of number- density evolution for the Early-type group. 



In lAbraham et al.l (120071 ) it is shown that evolution in the fraction of the stellar mass 
locked in massive early-type galaxies is produced in the interval 0.7 < z < 1.7. A model in 
which 0* for the Elliptical galaxies is constant to z ~ 0.6 and then evolve as 0* oc (0.4 + z)~ 2 
for higher redshifts also produces a good fit to the optical and NIR counts (see Figs. [231 
EH [251 and [26]). In this model the population of red elliptical galaxies has doubled since 
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Fig. 22. — Color evolution with redshift of the NIR colors for an Elliptical and an early 
spiral (with the parameterizations given in Tab. [8]). The two upper lines correspond to the 
J-Ks color whereas the other two correspond to the J-H color. The horizontal dotted lines 
correspond to the J-H~0.97 and J-Ks~1.84 colors of the fitted point for the slope change. 
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Table 9. Schechter Parameters for the Luminosity Functions 



Galaxy type 


M* AB {r') b 


4>* x 10~ 3 Mpc~ 3 


a 


E/SO 


-21.53 


1.61 


-0.83 


early Sp 


-21.08 


3.26 


-1.15 


late Sp 


-21.08 


1.48 


-0.71 


Im 


-20.78 


0.37 


-1.90 



a Considering Hq = 70 

b The characteristic galaxy luminosity given in the 
Sloan r' band in AB system. 




15 20 25 

B(Vega) 



Fig. 23. — Galaxy number count in B band taken from the literature. The lines correspond 
to the number counts predictions from the models described in the text. 
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J(Vega) 



Fig. 24. — Galaxy number counts in the J filter compared with data from other surveys. 
The lines correspond to number counts models described in the text. 
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H(Vega) 



Fig. 25. — Galaxy number counts in the H filter compared with data from other surveys and 
the number count models described in the text. 
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Fig. 26. — Galaxy number counts in the Ks filter compared with data from the literature 
and the galaxy counts models described in the text. 



z = 1, in good a greement w i th the increase of a factor of ~ 2 in the number evolution of red 
galaxies given in iBell et al.l (120041 ) . 



Finally, we have implemented a simple recipe to simulate downsizing in the elliptical 
population by maintaining 0* constant with redshift for the LF of bright galaxies. The 
results are compatible with our NIR galaxy counts and B-band counts from the literature in 
the case that <fi* is constant with redshift for red-ellipticals brighter than M* — 0.7 (~ —22.0 
in the Sloan r' band in AB system), decreasing the number densities for the bulk of the 
ellipticals as 0* oc (1 + z)~ 2 . 



6. Color analysis 

More information about the evolution of the galaxy populations could be obtained from 
color histograms. The separate number counts in each band at the magnitude ranges that we 
are sampling are less sensitive to the formation redshift (for values zf >= 4) or the e- folding 
timescale of the star formation than color histograms. Figure [28] shows the color-magnitude 
diagram builded with the ALHAMBRA data through the filter centered at 6130 A (F613) 
and Ks. The modelled evolutionary tracks for the 4 galaxy spectra considered in Tab. [8] are 
also displayed. Models with no evolution (top panel) and passive evolution (bottom panel 
) have been considered. As can be appreciated the evolved spectra produce better match 
to the data than the no evolved version, principally at the faint blue end which is better 
described by models considering passive evolution in the late Sp and Irr spectra. 

In Fig. [2H]it is shown the F613-Ks color histogram for different Ks magnitude bins. We 
have used an e- folding timescale r = 0.7 Gyr to describe the Elliptical galaxies in those 
plots, this longer timescale produce better fits to the red end of the color histograms than 
an instantaneous star forming event which over-predict the number of red galaxies. Both 
models produce similar results when fitting the galaxy counts in individual optical and 
NIR bands. The simulated histograms correspond to models where the population of E/S0 
galaxies decrease with redshift, the population of spirals stay constant, and the late type 
galaxies increase as as 0* oc (1 + z) 2 conserving the luminosity densities. The population 
of early spiral galaxies have been divided in two classes: one remain as in Tab. [8] whereas 
the amount of extinction have been doubled for the other. This try to avoid the fact that 
the discretization of the actual galaxy population in four classes tend to produce sharp 
histograms. 

As can be appreciated the models reproduce the overall shape of the data for bright 
Ks magnitudes. Nevertheless, at faint Ks magnitudes the models predict higher values in 
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14 15 16 17 18 19 20 21 

Ks (Vega) 



Fig. 28. — Color-magnitude diagram for the Alhambra data in the filter with effective wave- 
length 6130A (F613) and Ks. The tracks for the galaxy classes in Tab. [8] formed at zf=4 are 
also shown, for the passive evolution scenario top panel), and the no-evolution case bottom 
panel). Each galaxy track is computed for the characteristic luminosity given in Tab. [9j The 
solid tracks from top to bottom correspond to a simulated E/SO galaxy with single stellar 
population and e-folding timescales r = 0.5, lGyr. The dashed lines correspond to the early 
and late type spiral galaxies, and the dotted line to the Irr galaxy. The vertical dotted line 
and the diagonal one are the 5a detection levels. The two dashed lines join the spectra 
points for z=0.5 and z=1.0. The median color for each Ks bin is marked with crosses. 
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Fig. 29. — Color F613-Ks histogram (normalized for 1 square degree) for sources in different 
Ks bins. The histograms produced by the models have been convolved with a Gaussian 
kernel a = 0.2 in order to reproduce the photometric errors. The histograms computed for 
each population are shown in red for the E/S0, green for the early type spiral galaxies, and 
in blue for the late type spirals and Irr galaxies. 
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the color range 3 < F613 — Ks < 5. As could be inferred from Fig. [29] to obtain a better 
match to the data the number densities of early spiral galaxies formed in a shorter time- 
scale has to decrease with z, such fading of the spirals will lead to an under-prediction of 
the blue-band number counts unless that the number of star-forming increases at higher 
redshift. In Fig. [30] the F613-Ks histograms for the Ks bins: [16.5,17.5], [17.5,18.5] show a 
better concordance with the observed data. Those histograms correspond to a models where 
the number densities for the early spirals decrease with z as 0* oc (1 + z) -1 , the late type 
spirals number density remain constant with redshift, and number densities of Irr galaxies 
increase oc (1 + z) 3 . The luminosity density in not conserved within any galaxy class. Simila r 
results could be obtained using (f>* oc exp [(—Q/0)((1 + z)" 13 — 1)] (IBroadhurst et al.lll992l ). 
with Q=-l, and Q=3 to describe the number evolution of early spiral and Irr galaxies. This 
number evolution formulation avoid a step increase of Irr galaxies at high redshift. However 
in Fig. [301 a shortage of red galaxies is seen at about F613-Ks~ 5, this could be due to the 
fact that we only use a discrete number of galaxy parameterizations, for example is well 
know the exi stence of dusty sta rburst with the same red colors than the passive extremely 
red objects ( iDaddi et al.l 120001 ) . Covering a wider range in galaxy internal extinction or 
formation timescale will tend to smear the bi-modality present in the simulated histograms. 
As could be seen in Fig. [3H the number counts produced by this model in B + NIR filters 
also produce good fits to the observed data points. In this models the number counts at 
faint magnitudes will be dominated by the star-forming galaxies. The number count slope 
at faint magnitudes will be — 0.4(a + 1) (IBershadyl 120031 ) . being a the slope of the dominant 
luminosity function for M « M*. With this parameterization the slope of the number 
counts tend to 0.36 at the fainter end. 
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7. Summary 

We have presented galaxy counts in the J,H, and Ks filters covering an area of 0.45 
square degrees and an average 50% detection efficiency depth of J~ 22.4, H~ 21.3 and 
Ks~ 20.0 (Vega system). The depth reached, and the precision of the counts over a range 
of five magnitudes makes the data valuable for examining the change of the count slope 
reported in the Ks filter, and to extend this examination to the J and H bands. We find that 
a change in slope occur in each of the NIR bands in the range J=[18.5,19.5], H=[18. 0,19.0] 
and Ks=[l 7.0,18.0]. The NIR colors where the break in the galaxy counts are found imply 
that this change is related to the population of red galaxies at z~l. 

We have compared our number counts results with predictions from a wide range of 
number count models, concluding that in order to reproduce the described changes in the 
NIR slopes, a decrease in bulk of the population of red elliptical galaxies is needed. Good 
fits to the B-band and NIR counts are obtained with a parameterization for the number 
evolution of the elliptical population as <fi* oc (1 + z)~ 2 with no accompanying evolution in 
M*, corresponding to evolution in which the majority of ellipticals formed in spiral-spiral 
mergers. 

Performing a color analysis show that also the population of early spirals has to decrease 
at higher redshift i n order to describe the color distribution in r-Ks. Models using the 



parameterization of iBroadhurst et al.l (Il992l ) 0* oc exp [(—Q//3)((l + z)~$ — 1)], with Q=-3, 
Q=-l, and Q=3 to describe the number evolution of ellipticals, early spiral and Irr galaxies, 
with no number density evolution for late-spiral systems, produce good fits to the observed 
distribution, avoiding at the same time a high number of young systems at high z. A 
good match to the optical and NIR data is also obtained if the population of red-galaxies 
in the models remain constant to z ~ 0.6 and afterwards its number density decrease as 
<p* oc (0.4 + z)~ 2 , or if the number density of red-ellipticals is constant with redshift for 
galaxies brighter than M* — 0.7 (~ —22.0 in the Sloan r' band in AB system), decreasing as 
0* oc (1 + z)~ 2 for the bulk of the ellipticals. 

Alhambra is processing the data obtained in 20 medium-band optical and 3 NIR filters 
reaching high quality photometric redshift measurements (Az/(1 + z) < 0.03). Also an 
accurate classification by Spectral Energy distribution will be acquired. Those data will 
allow for the study of the evolution of the different galaxy types to 2 ~ 1, which will 
complement the results given in this article, disentangling what populations contribute to 
the number counts at different redshift inter vals. Also the stu dy of number counts for red 



galaxy populations, passive EROS or BzK (IDaddi et al.l 120041 ) galaxies will constrain the 



formation redshift and formation timescale for massive Elliptical galaxies. 
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Fig. 30. — Color F613-Ks histogram for sources in different Ks bins. The histograms have 
been normalized to f square degree. The modelled histograms have been computed for 
models in which the number densities for the early spirals decrease with redshift as <fr* oc 
(1 + z)^ 1 , the late type spirals number density remain constant with redshift, and number 
densities of Irr galaxies increase oc (1 + z) 3 . The number of E/SO galaxies decrease with z as 
specified in the labels. The histograms obtained with this models have been convolved with 
a Gaussian kernel a = 0.2 in order to reproduce the photometric errors. 
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Fig. 31. — Differential number counts in the B and NIR bands compared with a models in 
which the number densities for the early spirals decrease with redshift as 0* oc (1 + z)" 1 , the 
late type spirals number density remain constant with redshift, and number densities of Irr 
galaxies increase oc (1 + z) 3 . 
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A. Individual image combination 



As mentioned in the text, to combine the processed images we used SWARP (IBertin et al 



20021 ). With this software, individual images were projected into subsections of the final 
frame using the inverse mapping, in which each output pixel center was associated to a po- 
sition in the input image at which it is interpolated. With this schema, the code corrected 
at the same time the geometrical distortions in the individual images using the astrometric 
information stored in the image headers. 
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A.l. Estimating the relative transparency 

Using a filtered version of the SExtractor catalogs computed over the sky subtracted 
images, an accurate estimate of the relative transparency was computed by tracing the high 
S/N objects in all the images. The relative transparency values were used inside SWARP 
to scale the individual images to the same flux level in order to uniformize the zero points 
in the outer dither areas, this allowed to use 2MASS catalogs to calibrate the ALHAMBRA 
NIR photometry in the final images. 



A. 2. Astrometry calibration 

When computing the resampled version of the individual frames, SWARP uses the WCS 
information stored in the headers. In order to obtain a better matching of the individual im- 
ages, firstly the pipeline calculated the external astrometrical solution for a reference image. 
The individual images were then calibrated internally with respect to it, thus obtaining the 
equatorial coordinates from the reference image. In this paper the individual image with bet- 
ter transparency in a given pointing was used as reference. However, to get a better internal 
astrometry between the different filters, after completing a pointing in the 23 ALHAMBRA 
filters, the images with better FWHM and transparency in a set of selected optical filters, 
are combined to produce a deep image that will be used afterwards as reference. 



We have determined that the USNO-B1.0 catalog (IMonet et al.ll2003l ) provides an ad- 
equate number of objects to perform a high quality external astrometric calibration. We 
used our own code to match the sources with brighter apparent magnitudes in the reference 
image with those in the USNO-B1.0. Once a meaningful number of pairs was identified, the 
CCMAP IRAF task was used in two iterations to acquire the required astrometric solution 
with a 2nd order polynomial. A histogram of the external astrometric solution rms and the 
number of objects used for the final images of ALH08 is shown in Fig. [32J The median 
external astrometry rms is 0.12" ± 0.01 in RA and 0.11" ± 0.01 DEC. 

Having calibrated a reference image, the rest of the individual images were calibrated 
internally. The median internal rms in the astrometric solution for the OMEGA2000 data 
used in this paper is 0.06" in RA and DEC using a median of 160 objects, as shown in Fig[ 
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Fig. 32. — Histograms of the astrometry solution rms for the ALH08 images. (Top panel) 
External USNO-B1.0. (Bottom panel) internal. The small panels show the histograms of 
the number of objects that enter into the final astrometry fitting 



A. 3. Image co-adding 



Swarp allows the user to choose among several interpolation functions for inverse map- 
ping. For selecting the more appropriate kernel we analyzed the resulting final image FWHM 
and its pixel-to-pixel correlation. Tab. [10] shows the correlation values, between adjacent 
pixels and for pixel pairs separated by 2 pixels, in the final images obtained using different 
available interpolation functions. As can be seen in the table, the bilinear function produce 
a higher correlation which translate into an underestimation of the flux errors. Using the 
Lanczos-3 function the FWHM of the final image was improved by ~ 0.05" compared with 
the nearest neighbor interpolation, whereas the auto-correlation at 1 pixel remain acceptable 
~ 0.16 (when the images are combined using the average). The Lanczos-4 function did not 
decrease substantially nor the FWHM neither the correlation at 1 pixel, on the contrary it 
produced large artifacts at the bad pixels and image borders, so finally we have decided to 
take the Lanczos-3 function. 



-61 - 



Table 10. Values for the correlation between adjacent pixels when co- adding the images 

using different interpolation functions 



Nearest Bilinear Lanczos2 Lanczos3 Lanczos4 

Combining with median 

lpix +0.017 +0.186 +0.112 +0.074 +0.068 
2pix -0.035 -0.013 -0.072 -0.059 -0.101 



Combining with average 
lpix +0.035 +0.274 +0.177 +0.156 +0.130 
2pix -0.030 -0.018 -0.069 -0.110 -0.083 
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